arXiv:l 507.08097V 1 [math-ph] 29 Jul 2015 


Exact Solutions of the Heisenberg 
Ferromagnetic Equation 

F. Demontis* S. Lombardoj 
M. Sommacalj C. van der Mee* and F. Vargiu* 


Abstract 

In this article we construct an explicit (multi)-soliton solution for¬ 
mula for the Heisenberg ferromagnetic chain equation by using the 
Inverse Scattering Transform and the matrix triplet method. 


1 Introduction 


In this article we study the Heisenberg Ferromagnetic equation (HF), i.e., 


rrit = m A 


( 1 . 1 ) 


where m{z,t) G is a vector function satisfying m{z,t) —)■ 63 as 2 ; —)■ ±cx3 
and ||m(^,f)ll = ll^sil = 1 with 63 = (0,0,1)^, and A denotes the vector 
product. Here 2 ; denotes position and t time. 


Equation (1.1) describes the dynamics of the magnetization vector m 


of an isotropic ferromagnetic chain at the nanoscale in the absence of an 
external magnetic held. 

We underline that this article is the preliminary version of a much longer 
paper currently in preparation, which will account for all the proofs of the 
outcomes reported herein, and discuss in a more detailed and extensive way 


the exact solutions arising from our main result, i.e. formula (3.10), along 
with their interactions and classihcation. 


*Dip. Matematica, Universita di Cagliari, Viale Merello 92, 09121 Cagliari, Italy. 
iDept. of Mathematics and Information Sciences, Northumbria University, Newcastle, 
UK. 
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It is well known that (1.1) is integrable in the sense that it is possible 
to hnd a Lax pair associated with it (see [Ej). For later (notational) conve¬ 
nience, let us briefly recall that if 1/ is a 2 x 2 invertible matrix depending 
on position 2; G M, time t G M, and a spectral parameter A, the Lax pair 


associated to (1.1) is given by: 


V^ = AV= liA(m ■ (t)]V 

Vt = BV = ■ (j) — iX{T ■ (t)]V 


( 1 . 2 ) 


In (|1.2|) T is dehned as r = m A and cr is the column vector with entries 

0 -i'' 


the Pauli matrices 

'0 1 


CTi = 


1 0 


(72 = 


0 


0-3 = 


1 0 
0 -1 


Of course, the knowledge of the Lax pairs for (1.1) assures that the Inverse 


Scattering Tranform (1ST) (see [Il[21[T3]) can be applied to solve the initial- 
value problem for this equation, i.e.. 


rrit = m A rrizz, 
m[z,0) known 


(1.3) 


The hrst authors who applied the 1ST to the HF equation were Takhtajan 
and Zakharov [HIITj. 

The aim of this paper is twofold. The hrst goal is to present a (somehow 
novel) rigorous theory for the inverse scattering transform for the HF equa¬ 
tion. In particular, the direct scattering problem is proved to be well posed 
for potentials such that: 

1 . m{z) ■ cr has an almost everywhere existing derivative with respect to 
z with entries in 


2. m^{z) > — 1 for each 2; G M. 

Such conditions are less restrictive with respect to the usual Schwartz class 
hypotheses. For potentials satisfying 1. and 2. we establish the analyticity 
properties of eigenfunctions and scattering data. In order to obtain these 


results we derive a convenient set of dost solutions (see (2.10) below) which 


allows to study their asymptotic behavior at large A. The inverse scattering 
problem is formulated in terms of the Marchenko integral equations. These 
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are obtained by using an appropriate triangular representation of the Jost 
solutions (see (2.13) in Section]^. 

The second objective of this paper is to hnd an explicit multi-soliton 
solution formula for (1.1), which allows a classihcation of all its localized so¬ 
lutions, along with a description of their interaction. In order to achieve this 
goal, we will apply the matrix triplet method, which has proved successful 
in solving exactly - in the reflectionless case - several integrable equations 
[SllTlEllSlIlo]. The idea of this method is to represent the Marchenko kernel 
as (where {A, B,C) is a suitable matrix triplet) in such a way 

that the Marchenko integral equation can be solved explicitly via separation 
of variables. The solutions obtained in this way will contain nothing more 
complicated than matrix exponentials and solutions of Lyapunov or Sylvester 
matrix equations mu, hence can potentially be ’’unzipped” into lengthy 
expressions containing elementary functions. In the present preliminary pa¬ 
per we will show how to recover the famous HF one-soliton solution from 

( ra . 

The paper is organized as follows. In Section]^ we study the analyticity of 
the Jost solutions and scattering data, and determine their time evolution. 
Furthermore, we formulate the inverse scattering problem in terms of the 
Marchenko integral equation. In Sectioncombining the 1ST and the matrix 
triplet method, we get an explicit solution formula for (1.1). As a hrst 
example, we derive the explicit expression of the one-soliton solution. 


2 Direct Scattering Theory 


In this section we study the direct and inverse scattering theory associated to 


the hrst of equation (1.2). In particular, we analyze the analytic properties 


and the asymptotic behavior at large A of the Jost solutions and the scattering 
data, and formulate the inverse scattering problem in terms of the Marchenko 
integral equations. 

The proofs of the main theorems given in this section, i.e. Propositions 


2.3, 2.4 and Theorems 2.5, 2.11, can be found in [T6l ITT] . 


2.1 Jost solutions 

Let us dehne the Jost matrices J/(z, A) and <F(z, A) as those solutions of the 
linear eigenvalue problems = AT and = AT (A is the matrix dehned 
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in ( 1 . 2 )) satisfying the asymptotic conditions: 


A) = ('0(2;, A) •0(2;, A)) = + o(l)], 2 : -)■+oo, ( 2 . 1 a) 

$(z,A) = (0(z,A) 0(z,A)) =0^""3[J2 + o(1)], z^-cx). (2.1b) 


The columns 0 ( 2 ;, A), 0(^, A), 0 ( 2 ;, A), and 0 ( 2 ;, A) are called Jost functions. 
In the sequel, we also use the following notations: 


vI/(z,A) 


'rn^,x) 0“"(^,A)\ rnz,x)\ 

0""(z,A);’ ^ 1^0""(^,A) 


Then the differential equations = AT and 
be written as 


AT (cf. with ( 1 . 2 )) can 


=-iA(m ■ cr)T, (2.2a) 

=-iA(m ■ cr)T. (2.2b) 

It is then easily verihecQthat ^( 2 ;, A) and $( 2 ;, A) belong to the group SU{2). 
As a result. 


'J'ii(^)-^)* — ^ 22(^5 A), Ti 2 (^, A)* — —T 2 i(^, A), (2.3a) 

‘^’ 11 ( 2 ^ 5 '^)* = *^* 22(^5 A), $ 12 ( 2 ;, A)* = —$ 21(^5 A). (2.3b) 

Since the two Jost matrices are both solutions to the same hrst order linear 
homogeneous differential system, there exists a so-called transition matrix 
T(A), depending on A and belonging to SU{2), such that 

T(z, A) = T(z, A)T(A), A e M. (2.4) 


For A G M, we have 

fa{X) -b{X)\ 

{b{xy a{xyj^ 

where |a(A)0 -|- |&(A)0 = 1. We assume that a(A) y 0 for each A G M, i.e., 
we assume that no spectral singularities exist. 

In order to formulate the Riemann-Hilbert problem we need to establish 
the properties of analyticity as well as the asymptotic behavior of the Jost 

^Any square matrix solution 5(z) to the differential system = C{z)'E., where C{z) 
is skew-hermitian and has zero trace, has and detS independent of z G M. Here the 
dagger denotes the complex conjugate transpose. 
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solutions and the coefficients a(A) and 6 (A) at large A. To get these results, 
let us put mP = m — 63 . We can convert the differential systems (2.2) 
with corresponding asymptotic conditions (2.1) into the Volterra integral 
equations 

poo 

A) = -iX j dz ■ cr)T(A, A), (2.5a) 

J Z 

$(z,A) = +iA / ■cr)<h(A,A). (2.5b) 


As a result of Gronwall’s inequality (see Appendix of 0 ) we get for {z, A) G 
10)2 


11 ^( 2 ;, A)II < exp 1^1 A| ^ dA ||m°(A)||^ , (2.6a) 

||$( 2 ;, A)|| < exp |^|A| y dA ||m°(A)||^ , (2.6b) 

where we assume that m^{z) = m{z) — 63 has its entries in L^(]R). 

We can easily prove the following 

Proposition 2.1 If mP{z) = m{z) — 63 has its entries in L^(M), then the 
Faddeev functions e~''^^ip‘^"{z, A), A), and A) 

are analytic in X E C+ and continuous in X E C+ UM, while the Faddeev /unc- 
tions e^^^'fi''^{z,X),e^^^f)'^"{z,X),e~^^^4f^{z,X) and e~'^^^(j)‘‘"{z, X) are analytic 
in X E C“ and continuous in X E C“ U M. 

Proof. We give the proof only for the the Faddeev functions e~^^^'ip'"^{z, A), 
e-*^^'i/'d"(^, A), A), “( 2 ;, A) because the proof for the other Fad¬ 

deev functions is very similar. The Volterra integral equations can be written 
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in the form 


X) = i-i\ I dz A) 


0 / ;hA 


-iX dzm-{z)e-^^^ij‘^^{z,X), 

J Z 

POO 

X)=iX / dz A) 

J Z 

POO 

-iX dAe2*^("-")m+(A)e-*^V(A,A), 

J Z 

POO 

X) = -iX / dz A) 

J Z 

POO 


X) = l + iX I dz A) 


0 / jhA ^iXz^ 


— iX dzm^{z)e'‘^^iir^{z^X). 


(2.7a) 


(2.7b) 


(2.7c) 


(2.7d) 




Using Gronwall’s inequality, uniformly in (A, z) for A G C and z > zo> —oo, 


we obtain continuity in A G and analyticity in 


□ 


Taking the limits of (2.7a) and (2.7d) as z ^ —oo we get 


a(A) = l-iX dz [m°(A)e-*^V(5, A)+m-(A)e-*^>""(A, A)] , (2.8a) 

J —OO 

/ OO 

dz m°(A)e*^^'0'^"(A, A) — m'''(A)e*'^^'0“’"(A, A) , (2.8b) 

•OO ^ 

Thus a(A) is continuous in A G C+, is analytic in A G C+, and satishes 
a(0) = 1. In the same way we get 

/ OO 

dze^^^^ [m°(A)e-'^V“(5, A) - m+(A)e-'^"^/>“'’(A, A)] , (2.9a) 

■OO 

/ OO r n 

dAe-2*^" m°(A)e*^>“"(A,A)+m-(A)e*^^^""(A,A) . (2.9b) 

-OO ^ 

Thus 6(0) = 0 and 6 a (0) exists. 
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From (2.8) and (2.9) it is clear that no information is available on their 
asymptotics as A —)■ cxd. In order to get such information, let us derive a 
different set of Volterra integral equations. To do so we need to require that 

a. m{z) (7 has an almost everywhere existing derivative m'{z) ■ cr with 
respect to ^ which has its entries in L^(M) 

b. m^i^z) > — 1 for each z G M. 

Under the above hypothesis, after straightforward calculations we get 


D{z)'^{z,X) = e 


_ „jA2(T3 


A), 


( 2 . 10 ) 


where 


( 2 . 11 ) 

which is a matrix of determinant ^(1 + ms^z)). Since we have required that 
m^{z) > — 1 for each z G M this matrix D{z) is invertible and its inver s^ is 
bounded in z G M. We may therefore apply Gronwall’s inequality to (2.10) 
and hnd that 


il^(uA)||< 


^/TTm^ 


exp 


dz ||(m'(A) ■ cr) 


Remark 2.2 In the same way and under the assumptions a and b, adapting 
the procedure above presented to the lost matrix <F(^, A), we get 


D{z)^{z,\) = e 


— „*Azo-3 


+ 


^;ggU(2-2)a3^/(^)^(^,A). (2.12) 


We may therefore apply Gronwall’s ineguality to (|2.12|) to obtain 

1 


II®(UA)||< 


^yi + msiz) 


exp 


dz \\{m'{z) ■ 


^Under the first condition, m{z) is absolutely continuous in z G M. Hence its pointwise 
values make sense and hence it makes mathematical sense to assume that, in addition, 
7713(2) > — 1 for each 2 G K. 

^That is, D{z)~^ = ^ niatrix of norm (1 + 7773)“^/^. Note that 

D{z) —>• U as 2 —>■ ±00. 
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Equations (2.10) and (2.12) allow us to prove that the analyticity and con¬ 
tinuity properties of the dost solutions extend to the closed upper and lower 
half-planes. In other words, the dost solutions and the coefficient a(A) have a 
hnite limit as A —)■ cxd from within the closure of its half-plane of analyticity, 
while 6(A) vanishes as A —)■ ±cx). In order to prove these results we need to 
hnd a “suitable” triangular representation for the dost solutions. We have 
the following: 


Proposition 2.3 There exist an auxiliary matrix function K{x,y) such that 


'^{z, A) = + / dz K{z, z)e 


iXza^ 


(2.13) 


where H(z) is a matrix function satisfying H(z) = a 2 H(z )*a 2 , H{z) —)■ I 2 
as z ^ -|-oo, dz \\K{z, A)|| converges uniformly in z 


can be found in [151 n]. Analogously we have 

the following 


The proof of Proposition 2.3 


Proposition 2.4 There exists an auxiliary matrix function N{z,y) such 
that ^ 

^{z,X) = H{z)e^^^^^+ f dzN{z,z)e^^^^^, (2.14) 

J —00 

where H{z) is a matrix function satisfying H{z) = a 2 H[z)*a 2 , H {z) ^ h 
as z ^ — 00 , and Jl^dz || A/'(2;,;S)|| converges uniformly m z G M. 

We remark that the auxiliary matrices K{z, A) and N{z, A) have to satisfy 
the symmetry relations 

N,{ 

-N 2 

(2.15) 

where Ki{z,y), K 2 {z,y), Ni{z,y) and N 2 {z,y) are scalar functions. 
Finally we can prove (see [151 HU for the proof) the following: 

Theorem 2.5 7/1) m^{z) = m{z) — 63 has its entries in 2) m{z) a 

has an almost everywhere existing derivative m\z) ■ cr with respect to z 
which has its entries in and 3) that m^{z) > — 1 for each 2 ; G M, 

then for each z & ^ the functions A), A), A) 

and e^^^(f)'^"'{z, A) which are analytic in \ E C"*" and continuous in \ E 


z,y)* N 2 {z,y)\ 

{z,yy Nyz,y)J^ 


K{z,y) = 


Ki{z,y)* K 2 {z,y) 

-K2iz,yr Kyz,y) 


N{z,y) 





U M (see Proposition 2.1) have a finite limit as A —)■ cx) from within 
the closure of C+. Analogously, for each 2 ; G M the Faddeev functions 
A), A), A) and e~^^^(j/"'{z, A) which are analytic 


in X G C and continuous in X G C UM (see Proposition 2.1) admit a finite 
limit as A —)■ cxD from within the closure ofC~. Moreover, the coefficient a(A) 
has a finite limit when A —>■ 00 from within C"*" while 6 (A) may not admit 
analytical continuation outside the real line and 6 (A) —)■ 0 when X —)■ ±cx). 


2.2 Scattering Data 

In this subsection, for the sake of completeness, we introduce the scattering 
matrix and the scattering coefficients. 

From now on, we assume that the coefficients a(A) introduced in the pre¬ 
ceding section is such that a(A) 7 ^ 0 for each A G M, i.e., there are no spectral 


singularities. We can write the identity (2.4) as the following Riemann- 
Hilbert problems: 


z,X) '(^{z,X)) = {f){z,X) 4>{z,X)) 



(fi(z, A) (/)(z, A)) = ((/)(z, A) fi(z, A)) 


,a(A)* 

Putting F_{z, A) = [(j){z, A) fi{z, A)) and F^{z, A) = 
obtain the Riemann-Hilbert problem 

F.{z,X) = F+{z,X)a 3 S{X)a^, 

where the scattering matrix S{X) is dehned by 

_ fnx) Rix)\ 
T{X))- 


b(A)\ 

a(A) 

(2.16a) 

HP) J 

b(A) \ 

■F ■ 

(2.16b) 

xixrJ 

{'fi{z,X) 

0(z, A)), we 


(2.17) 


We write 


T(A) = 


R{X) = 


6 (A) 


m = 


6 (A)^ 


(2.18) 


a(A)’ a(A)’ a(A) ’ 

where T, R, and L are called the transmission coefficient, the reflection coef- 
hcient from the right, and the reflection coefficient from the left, respectively. 


Equations (2.16) then imply that 

R(A)t = a3R(A)-V3, 


A G 
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Thus S'(A) is (Ta-unitary and has determinant a(A)*/a(A). Also, S'(A) — 
e “*"/2 as A —)■ ±cxD for a suitable complex number of modulus 1. We 
easily derive the Fourier representations 


F+{z,X)e 


-iXzas _ (Hi{z) H2 {z) 

H2{z) H,{z) 


+ 


ds e* 


iXs (Fi{z,zZr s) N 2 {z,z-s) 


F_{z,\)e 


K2{z,z + s) Ni{z,z-s)J' 

-iXzaa _ I ~H2 {.z] 

H2{zy H^zy 

-iXs ( Ni{z,z-sy -K2{z,z + sy 
-N2 {z,z - s)* Ki{z,z + sy 


ds e~ 


where 


(2.19a) 


(2.19b) 


ds {\Ki{z, + s)| + \K 2 {z, + s)| + \Ni{z, 2 ; - s)| + \N 2 {z, 2 ; - s)|) 


converges uniformly in z G M. 


The scattering data associated with the hrst of equation (1.2) are: 


1 . one of the reflection coefficients, 

2 . the poles of the transmission coefficient T(A) (or of T(A*)*). We call 
such poles the discrete eigenvalues in the upper half-plane C"*" (or in the 
lower half-plane) and denote them by iaj (or by —ia*) for j = 1,..., N. 


3. a set of constants Nj (Nj) for j = 1,..., associated to the discrete 
eigenvalues iaj {—ia*) j = 1,N in the upper half-plane (lower half¬ 
plane). These constants are called the norming constants. 


By using elementary arguments of complex analysis it is possible to prove that 
if there are no spectral singularities then the number of discrete eigenvalues is 
finite. We present, for the sake of simplicity, the scattering theory assuming 
that each pole of the transmission coefficient has multiplicity equal to one 
(the general case can be treated as shown in [5]). The way to construct 
the norming constants is standard (see HEIE]). However, for the sake 
of completeness, we prefer to insert their construction. Let us assume that 
there are finitely many poles iai ,..., ioAr of the transmission coefficient T(A) 
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in the upper half-plane C"*", all of which are assumed to be simple. We let 
stand for the residue of T(A) at A = i.e., 


Ts = lim (A - m,)T(A) = lim —^ , 
x^tas a[\) — a[ias) a[ias) 

where the dot indicates differentiation with respect to A. We then introduce 
the norming constants Ng by 


Ts4>{z,ias) = iNs'ijj{z,ias), 


s = l,2,...,N. 


(2.20a) 


By the same token, T(A*)* has the simple poles —ial ,..., —ia*j^ in C , all of 
them simple. The corresponding norming constants Ng are dehned by 


T* (j){z, -ial) = -iNgi){z, -ial), s = 1,2,..., W 


(2.20b) 


The next proposition shows how the norming constants introduced in the 
upper half-plane are related to those dehned in the lower half-plane. 

Proposition 2.6 The norming constants satisfy the following relations: 

Ng = -{Ng)* 


Proof. By applying the triangular representations to (2.20a) and (2.20b) 
we get the pair of equalities 


H2{z)\ r f N2 {z, z-w) 

Jo [n,{z,z-w) 


= tNg 


Jo [K 2 {z,Z + w) 




= -iN 


-H2iz''*^ 
HAz) 


Ni{z, z — w)* 

-N 2 {z, Z - w)" 

K 2 {z, Z + wf 


+ dwe . 

Jo \Ki{z,z + w) 


Taking the complex conjugate of the hrst equation and premultiplying the 
result by (o), we obtain the second equation, provided Ng = —{Ng)*. □ 
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2.3 Marchenko Equations 

In this subsection we formulate the Marchenko integral equations and estab¬ 
lish the connection between the solutions of these equations and the solution 
of the HF equation. We refer the reader to [ini HB E] for the details on the 
derivation of (2.24) below. 

In order to derive the Marchenko equations we need the following 


Proposition 2.7 Suppose that hypotheses a and b in Subsection \2.1\ hold, 
and suppose there are no spectral singularities. Then there exist functions p 
and I in L^(M) such that 


R{X)= ! L{X)= I dwR^'^l{w). 


iXw 7 


( 2 . 21 ) 


We refer to [iniin] for the proof. 
We have the following 


Theorem 2.8 The auxiliary function K{x,y) which appears in (2.13) sat¬ 
isfies the following integral Marchenko equations 


K{x,y) + H{z)n{zX-y) + 
where 
fl(w) = 


duK{z,u)fl{u + y) = 02 x 2 , ( 2 . 22 ) 


0 uj{w) 


—ui{wy 


0 


j , with uj{w) = p{w) e , 

' et 


(2.23) 


and p{w) is the Fourier transform of the reflection coefficient (see (2.21)J. 

Recall that H(z) G SU(^2). Writing K(z,y) = H(z)L(z,y) we can 
convert (2.22) into the traditional Marchenko integral equation. By following 
the same proof as in the focusing AKNS case EE], we hnd that the integral 
equation 


L{z,y) + fl{z + y) + / du L{z,u)fl(u + y) = 02x2, 


(2.24) 


is uniquely solvable on the space L^{z, -l-oo)^^^. We observe that analogous 
Marchenko equations are satished by the auxiliary function N{x,y) which 
appear in (2.14). More precisely, we have the following 
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Theorem 2.9 The auxiliary function N{x,y) which appears in (2.14) sat¬ 
isfies the following integral Marchenko eguations 

N{x,y) + H{z)Cl{z+ y) + j du N{z,u)tt{u + y) = (2.25) 

J —OO 

where 

Q(w) = 


0 

—u)(w) 


u{w) 

0 


con ui{w) = l{w) + ^2 ^ , (2.26) 


S=1 


and l{w) is the Fourier transform of the reflection coefficient (see ( 2.21 [ ) ). 

Setting N{z,y) = H{z)P{z,y), we can convert (2.25) in the following inte¬ 
gral Marchenko equations 


P{z, y) -h Ft{z + y)+ j duP{z, y)Q{u + y) = 62x2 • 


(2.27) 


From now on we get our results working with the equations (2.22) and 


(2.24). We remark that similar results can be established starting with equa¬ 
tions (2.25) and (2.27). 

For later convenience let us introduce the following notations 


K(z) = 


dyK{z,y), L{z)= dyL{z,y), 


(2.28) 


where K{z,y) and L{z,y) satisfy the Marchenko integral equations (2.22) 
and (2.24), respectively. Using the Volterra equation (2.5a) and the asymp¬ 
totic relation \I/(z, —)■ J 2 as z —>■ -foo, we get from the triangular 

(2.29) 


representation (2.13) 


h = ^(^, 0) = H{z) + K{z) = H{z) 


I2 + L(z) 


where L{z) = H-\z)k{z). 

The relationship between the Marchenko integral equation and the solu¬ 
tion of equation (1.1) is immediately clarihed by the following 


Proposition 2.10 Under the hypoteses a and b of Subsection 2.1, we have: 


m{z) ■ (T = H{z) a3 H(z) I2 + L{z) 

We refer the reader to [T6| ITT] for the proof. 


-1 


CTS 


I 2 + L{z) 


(2.30) 
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2.4 Time Evolution of the Scattering Data 


In this subsection we derive the time evolution of the scattering data. We 
shall arrive at the same time evolution as for the NLS equation. 

Recall the AKNS pair is given by (1.2). Suppose V(z, t; A) is a nonsingular 
2 x2 matrix function satisfying 


K = AI/, Id = BV, 


where V does not need to be one of the dost matrices. Then there exist 
invertible matrices Cq, and Cq,, depending on {t, A) but not on z, such that 
^ and d) = VC^\ Thus 

= VtC^^ - VC^^[Cq.]tC^^ = BVC^^ - VC^^[Cq,]tC^^ 

= Bm 


implying 

[Cq,]tC^^ = (2.31) 


Here the left-hand side does not depend on x, whereas the right-hand side 
only seemingly depends on x. We may therefore allow z to tend to -|-oo 
without losing the validity of (2.31). Since B ~ —2iA^(T3 and T ~ 
z —)■ -fcx), we obtain 


^iXza^ 


as 


[Cq,]tC^^ = -2zAV3. 


(2.32) 


The same result can be obtained for the other dost matrix <l>(z, A). Using 
that T = ^T, we get 


Tt = (<I)-^T)t = 

= $-i(RT - ^Cq,]tC^^) - <h-'(R$ - <I>[C'$]tC'^i)<h-iT 
= - T[Cq,]tC^^ - 

= 2zA^ {Tcr3-cr3T}, 


so that 

Hence, a(A) and T(A) do not depend on t, whereas 

R{X, t) = e-^*^'‘R(A, 0), L(A, t) = e^*^'T(A, 0). 


Differentiating (2.20a) with respect to t, we obtain 


Ts4>t{z,ias) = iNs'iptizRas) + i[Ns]t'ip{z,ias). 


(2.33) 

(2.34) 
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Using ( 2.31[ ) for $ and \I/ and using (2.32), we get 
Ts {B{ias)(l>{z,ias) - 

= iNs {B{ias)'il^{z,ias) + 2m^'0(2;, m^)} + i[Ns\t'ip{z,ias). 


Using (2.20a) again we obtain 


[iV,]t = -AtalN,. 


Taking into account that Ng = — [A^s]* (see Proposition 2.6), finally we get 
Ng{t) = Ns{t) = (2.35) 


2.5 Inverse Scattering Transform. 

Having presented the direct scattering problem (consisting of the construction 
of the scattering data when m{z, 0) is known), the inverse scattering problem 
(amounting to the construction of ra{z) when the scattering data are given), 
and the time evolution of the scattering data associated to the hrst of equation 


(1.2), we can discuss how the 1ST allows us to obtain the solution to the initial 


value problem for (1.1). 


Using the initial condition m(z, 0) as a potential in system (1.3), we 
develop the direct scattering theory as shown above and build the scattering 
data. Successively, let the initial scattering data evolve in time in agreement 


with equation (2.33)-(2.35). The solution of the Heisenberg equation is then 
obtained by solving the Marchenko equation (2.24) where the kernel fl{w) is 


replaced by Q,{w,t) (i.e., taking into account (2.33), (2.34), and (2.35)), and 


then using relation (2.30). 


Gauge transformation. In this paragraph we show how the gauge 
transformation (see miin]) between the solutions of the Heisenberg equation 
and the solutions of the Nonlinear Schrodinger equation is determined by 
the matrix H{z) {H{z)) introduced in the triangular representation (2.13) 


(CT). 

Let us denote by •^)) the dost matrices of the Zakharov- 

Shabat system, and by T(z, A), <I)( 2 :, A) the dost matrices of the scatter¬ 


ing problem (1.2) associated to the HF equation, both corresponding to 
the same scattering data. We have already established the following rela¬ 
tions; H{z) = [h + dwL{z,w)~\ ^ (see (2.28) and (2.29)) and H{z) = 


f Z 

[I 2 + dwN{z,w)]~^ (see ( [Tiil ), pTS] ), ( [T^ and p(^ ). 


Then we have the following: 
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Theorem 2.11 The solution of the initial value problem (1.3) are expressed 
in terms of the Jost solutions of the Zakharov-Shabat system as: 


m{z) ■ a = Q)(T 2 ,^zs{z, 0) = $^ 5 ( 2 ;, 0)a3^zs{z, 0). (2.36) 

The proof can be found in [161 E] ■ 


3 Matrix Triplet Method 

In this section we construct an explicit soliton solution formula for equation 
). We apply the same technique successfully used in [61171181191 HH] 
to solve the NLS, the cmKdV, the sine-Gordon, the discrete NLS, and the 
Hirota equations, respectively. Furthermore, we use the triplet method to 
get explicit expressions for the Jost solutions in the reflectionless case when 
the corresponding scattering data are specified. 



3.1 Explicit soliton solutions for equations (1.1). 

We want to restrict ourselves to the case R{\) = 0. In this case the expression 
for f2i(t(;;0) is given by (2.23), with p = 0. In particular, we can treat the 
situation where the discrete eigenvalues are not necessarily simple 0 by 
introducing the function 


N n-j-l 




s=l s=0 


W 

s! 


(3.1) 


rii — l 
0 


t=l 


where (3.1), {i are the discrete eigenvalues in C"*", |{A(js(t)}gi, 

are the corresponding norming constants, and Uj is the geometric multiplicity 
of i Oj. 

To recover the solution of (1.3) we follow the three steps indicated below: 
a. Assume that the discrete eigenvalues and the norming con¬ 


stants \ 


N 


i=i 


are given, with N denoting the number of 


discrete eigenvalues in C’*', and n,- denoting the geometric multiplicity 
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of ittj. By using the scattering data we introduce the matrix function 


' ’ ' ' Ef..ES 0 


(3.2) 


where the norming constants evolve in time according to (2.35). 


b. Solve the integral Marchenko equation (2.24): 


L{z,y;t) + fl{z + y;t) + / du L{z,u-,t)fl{u + y;t) = 02x2- 


where u > z, and the kernel fl{z,y) is given by ( |3.2[ ). 
c. Construct the potential m{z;t) by using formula([2.30[): 


m{z]t) ■ cr = I 2 + L{z]t) 


^ -1 




I 2 + L{z] t) 


where L{z) = J^dyL{z,y). 

An analogous procedure can be followed by using the Marchenko equation 


(2.27). 


Let us follow the above procedure, momentarily disregarding the time 
dependence. We will then show how to take into account also the time 
dependence. 

It is well known [T^l E] that it is possible to factorize a matrix function 


in the form (3.2) by using a suitable triplet of matrices. More precisely, 


suppose {A, B, C) is a triplet of matrices such that all the eigenvalues of the 
q X q matrix A have positive real parts, B is 2q x 2, and C is 2 x 2q. Let 
Q (w) he dehned as 


0 


Z^s=l Z^s=0 s! ^ 
0 


a(w) I ^ ^ 

2^s=l Z—is=0 '^'js s! ^ 


= Ce-'^^B. 
(3.3) 


Equation (3.3) can be written as (see (2.23)) 


f2(w) = 


0 

-oj{wy 


u:[w 
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where the scalar function lj{w) is dehned as u{w) = with A a. qx q 

matrix having only eigenvalues with positive real part, B a qxl matrix, and 
C a 1 X g matrix. Furthermore, let us assume that the triplet {A, B, C) is a 
minimal triplet in the sense that the matrix order of A is minimal among all 
triplets representing the same Marchenko kernel by means of (3.3) 

We then dehne 


Then 

where 


A 0 


qxq 


qxq 

A'f 


V = 


t2 _ / Ogxi ^ 

' -C^ %X 1 


C = 


c 

Olxq 


Olxg 

St 


dze-^^BCe-^^ = 




qxq 

-Q 0, 


N = 


dz e 


-zA 


BB^e-^^\ 


N 

qxq 


Q= I dze-^^^C^Ce 


.-zA 


(3.4) 


Then V is the unique solution of the Sylvester equation AV + VA = BC. It 
follows from the minimality of the triplet [A^B^C) [3l Sec. 4.1] that N and 
Q are positive Hermitian matrices. Thus V is invertible and 


-p-i _ / Ogxg Q 

' ' AT-i 


0 


qxq 


(3.5) 


Substituting the expression for the kernel (3.3) into equation (2.24), we 
arrive at the following Marchenko equation 


L{z,y)+Ce-^^+y^^B+ / duL{z,u)Ce-^^+y'^^B = Q 


2x2- 


(3.6) 


Equation (3.6) can be solved explicitly via separation of variables. In fact, 
looking for a solution of the following form 

L{z,y) = -F{z)e-y^B, 

after some straightforward calculations, we find 


L{z,y) = 


(3.7) 


provided the inverse matrix exists for each 2 ; G M. Actually, the inverse 
matrix exists thanks to the unique solvability of (3.6). Finally, in order to 
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reconstruct (1.3) we have to integrate (3.7) with respect to y, obtaining the 
explicit formula 




(3.8) 


By means of this latter expression we can write the right-hand side of (2.30) 
in explicit form, and recover the components mj{z) of the vector m{z). 

In order to introduce the time dependence we have to take into account 


the time evolution of the scattering data expressed by (2.33)-(2.35). Dehning 
the (reflectionless) Marchenko kernels as follows: 


u{w; t) = t) = (3.9) 

we may replace the matrix triplet (A, B, C) by (A, 5, in such a way 

that (2.33), (2.34) and (2.35) are satished {A contains the discrete eigenvalues 


which are time independent, and C the norming constants). 


As a consequence, the explicit right-hand side of (2.30) is given by: 


L{z-t) = -C{t)e-^^[hq + 


(3.10) 


where 


m = 


Jgxl 


- Oq 


C{t) = 


Olxq 


Olxg 

St 


V{t) = 


( 1 = [ dze-^^^ 

\ Qy'^) Ugxg/ Jo ^ ' 


We remark that formula (3.10) allows us to compute explicitly the functions 


mi(z,t) satisfying the HF equation (1.1). The following example illustrates 


this fact in the most elementary case: the one-soliton solution. 


Example 3.1 (One-soliton solution.) The easiest example is the one-soliton 
solution which is obtained by taking A = (a) with p = Re a > 0 , B = (l), 
and C = (c) with 0 7 ^ c G C. Then Q = (|cp/2p), N = (l/2p). Hence, 


by substituting (3.10) in the RHS of (2.30), after long but straightforward 
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calculations, we get 
Icl^ 

m^{z) = 1 + 


^-Spz 


c 

2 

e-4pz 

C 

2 

e-4pz 

a 

P l + (|c| 

2/4p2)g-4p2 

a 

P li + (W 



m+{z) 


1-2 
„ c 


c 

2 

^ 2pV|a| 

2 

a 

P [l + (|c| 

p/4p2)e“^P^]^ cosh^[2p(2 

- ^ 0 )] ’ 


1 


e- 4 pz 


—4pz 


^—2az 

1 

c 

2 

1 - 

1 

l + (|c 

2/4jg2)g-4p2 

2pa 1 J- ( c 

2/4p2)g-4p2 


C 

m_{z) = —2- 
a 

solution is obtained by replacing the matrix triplet (a, 1, 
in ( |3ll| ). 


(3.11) 

(3.12) 

(3.13) 

c) by (a, 1, 


3.2 Reconstruction of the Jost solutions. 

By using the same notations introduced in the subsection above, let us com¬ 
pute the Jost solution ^/(z, A) by substituting the solution of the Marchenko 
equation (2.22) (see formula (|3(T)) into (2.13). We get 


^{z, = H{z) - H{z)Ce-^^ [hq + / dy 

J Z 

POO 

= H{z) - H{z)C [V + / ds 

Jo 

= H{z) - H{z)C [V + W(A), 

where —MW(A) -|- iAW(A)cr 3 = BC has a unique solution]^ Since <73 and B 
anticommute, and A and <73 commute [see ( |3.4[ )], we have 

POO POO POO 

W(A)= / dse-^^Be^^^^^ = / = / ds 

Jo Jo Jo 

= {iXas + A)~^B = -iasiXhq - ia^Ay^B = i{Xl2q - ia^Ay^Bas. 


^Note that cr(yl) n {iA, —iX} = 0 for A S K. 
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Consequently, 




H{z) 

H{z) 


hq -iC[V + ■' (M 2 , - ^cT3A)-^&a3 

I 2 , - iCe-^^ [hq + e-^^iXh, - tasA^^Bas 


Taking the limit as 2 ; —>■ — 00 , we obtain 


a(A) 


lim e 

\z^—oo 


— 2iXz 


b(xy 


lim e‘^^H(Xy 

Z^ — OO 

a{xy 


h-iCV ^(Xhq- ia^A) 


provided p(w) is invertible. Equations (3.4) and (3.5) imply that 6(A) = 0 
[reflectionless case]. Moreover, we And the following: 


a(A) = 1 - iCQ-\XIq + iAy-^C\ (3.14a) 

a(A)* = 1 + iB^N-\XIq - iA)-^B. (3.14b) 


Taking the complex conjugate transpose, we get the alternative expressions 


Thus a(oo) 


a(X) = 1 - iB\XIq + iAy-^N-^B, 
a(Xy = 1 + iC(XIq - iA)-^Q-^Cy 


1. Taking determinants in any of (3.14), we getj^ 


_ det(AJg - iA) 

det(A/g + iA'^) 


(3.14c) 

(3.14d) 
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